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1 Introduction 



Since its discovery in 1975 [T], the instanton solution of the Yang-Mills equations has 
been widely studied, both in its mathematical properties and its phenomenological 
applications [21 [3l |31 13 13 13 E] . Many insights have been obtained in QCD through 
the use of instantons (see, e.g., the review [9] and references therein), and even if it 
is known that they cannot provide the framework for a complete understanding of 
strong interactions, it is interesting to investigate instanton effects as they contribute 
to the nonperturbative dynamics. 

Among the various open problems in QCD, soft high-energy hadron-hadron 
scattering is known to be of nonperturbative nature, and so it is worth studying what 
the consequences are of instantons on the scattering amplitudes. In the approach 
from the first principles of QCD [T^, such amplitudes are related to properties of 
the vacuum, namely the correlation functions of certain Wilson-line and Wilson- 
loop operators. In particular, in the case of meson-meson scattering, the scattering 
amplitudes can be reconstructed, after folding with the appropriate mesonic wave 
functions, from the scattering amplitude of two colour dipoles of fixed transverse 
size; the latter are obtained from the correlation function of two rectangular Wilson 
loops, describing (in the considered energy regime) the propagation of the colour 
dipoles 112 ^ [T3 l I14 ^ I15j. The physical, Minkowskian correlation functions can be 
reconstructed from the "corresponding" Euclidean correlation functions [161 [T71 [T8l 
\T9\ \20\ EH [22], and so it has then been possible to investigate the problem of soft 
high-energy scattering with some nonperturbative techniques available in Euclidean 
Quantum Field Theory [23l [Ml |25l [211 [271 [28] . 

In this paper we shall derive a quantitative prediction of instanton effects in 
the Euclidean loop-loop correlation function, relevant to the problem of soft high- 
energy scattering. This problem has already been addressed in Ref. [23], using the 
so-called Instanton Liquid Model (ILM) [8], but the result reported in that paper 
contains a severe divergence, apparently not noticed by the authors. In this paper 
we critically repeat the calculation, obtaining a well-defined analytic expression, 
which can be compared with the lattice data presented and discussed in Ref. [28j . 
In this paper we also calculate the instanton-induced dipole-dipole potential from 
the correlation function of two parallel Wilson loops [29j , and we compare the results 
with some preliminary data from the lattice. 

The detailed plan of this paper is the following. In Section [21 we briefly recall the 
main features of instantons in Yang-Mills theory, we shortly describe the relevant 
aspects of the ILM, and we discuss the method we will actually use in our calcula- 
tions, using as an example the expectation value of a rectangular Wilson loop and 
the instanton-induced gg-potential [H 130]. In Section [3l after a brief review of how 
high-energy scattering amplitudes can be reconstructed from the correlation func- 
tion of two rectangular Wilson loops at an angle 6 in Euclidean space, we evaluate 
this same correlation function and we critically compare the result with the calcu- 



1 



lation of Ref. [53]. We then compare our quantitative prediction with the lattice 
results of Ref. [28]. In Section [5] we calculate the correlation function of two parallel 
Wilson loops, from which we derive the instanton-induced dipole-dipole potential, 
that we compare with some preliminary data from the lattice. Finally, in Sectionj5] 
we draw our conclusions. Some technical details are discussed in the Appendices- 



2 Wilson loops in the ILM 

Instantons (anti-instantons) are self-dual (anti-self-dual) solutions with finite ac- 
tion of the classical Yang-Mills equations of motion in Euclidean space [1] . We will 
be interested only in the solutions with topological charge q = +1 and q = —1, 
which we will refer to as the instanton (I) and the anti-instanton (I), respectively. 
In the case of SU{2) Yang-Mills theory, the instanton solution in regular covariant 
(Lorenz) gauge reads, for the standard colour orientation, 

^>''^'^ = -lw=WT^^ ^. = ^^T' « = 1'2,3, (2.1) 

where and p are free parameters that correspond to the position of the center and 
the radius of the instanton, respectively, ija^u are the self-dual 't Hooft symbols [5], 
and a" are the Pauli matrices. The anti-instanton solution A1^(x;z,I) is obtained 
by replacing rjafiu rja^u in Eq- (12.11) . where the latter are the anti-self-dual 't 
Hooft symbols |5]. Applying a SU{2) transformation on Eq. ()2.ip . UA^U'^ 
with U G SU{2), we obtain another solution with a different colour orientation. 
Instantons in SU{Nc) Yang-Mills theory can be obtained by embedding the SU{2) 
solution (|2.ip in a SU{2) subgroup of SU{Nc), and it has been shown that these are 
the only solutions with \q\ = 1 [7]. We will take the standard-oriented solution to 
be the embedding of the (standard) SU (2) solution in the upper-left 2x2 corner of 
an Nc X Nc matrix; the other solutions are obtained by applying an SU{Nc) colour 
rotation on this one. 

Early attempts at a description of the QCD vacuum in terms of an ensemble of 
instantons and anti-instantons adopted the picture of a dilute gas [6], but suffered 
from severe infrared divergences due to the presence of large-size instantons. The 
proposal of [8j was instead that the pseudoparticles form a dilute liquid: the di- 
luteness of the medium allows for a meaningful description in terms of individual 
pseudoparticles, while the infrared problem is solved assuming that interactions 
stabilise the instanton radius. In particular, this model assumes that the vacuum 
can be described as a liquid of instantons and anti-instantons, with equal densities 

^In this paper we will deal only with the Euclidean theory, and so we will understand that the metric is 
Euclidean and that expectation values have to be taken with respect to the Euclidean functional integral, 
except where explicitly stated. 
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ni = nj = n/2, and with radius distributions strongly peaked around the same 
average value: in practice, one performs the calculations using a 5-like distribution, 
i.e., fixing the radii of the pseudoparticles to a given value p. Using the SVZ sum 
rules |3H I32j to relate vacuum properties with experimental quantities, the phe- 
nomenological values for the total density and the average radius are estimated to 
be [9] n ~ lfm~^, p ~ 1/3 fm, which lead to a small packing fraction nVi ~ 1/20 
[here Vi = {tt"^ /2)p'^ is the volume of a 3-sphere of radius p, which is approxi- 
mately the space-time volume occupied by a single pseudoparticle] ; this indicates 
that the instanton liquid is indeed fairly dilutqj. We refer the interested reader 
to [9] and references therein for a more detailed discussion of the model, both on 
its theoretical foundations |33l |3^ and on the development of numerical methods 
of investigation ^35j . In the rest of this paper we will make use only of the general 
picture described above to give estimates of the relevant quantities. 

To calculate the expectation value of an operator O one needs to average over 
the ensemble with some density function V = P({2:}, {[/}), where {z} denotes 
collectively the coordinates of the pseudoparticles and {[/} their orientations in 
colour space; in principle, such a function should be determined by properly taking 
instanton interactions into account. Starting with an ensemble of A^i instantons 
and Ni anti-instanton^ in a space-time volume V ^ we obtain the desired result 
by taking the "thermodynamic limit" Ni,Ni,V — )• oo, while keeping the densities 
m = Ni/V and nj = Ni/V fixed: 

(O)iLM = lim / d^"" zd''UV{{z},{U})0{{z},{U}). (2.2) 

rij J fixed 

Here ©({z}, {C/}) denotes the operator O evaluated in the field configuration gen- 
erated by the pseudoparticles, and dU is the normalised Haar measure, / dU = 1. 
Since the medium is dilute, the field configuration is approximately equal to the 
superposition of single (anti-) instanton fields, 

N 

A^{x) = Y,U{i)A^{x;z,,cji)U\i), (2.3) 

i=l 

where A'" = A^i -|- A'j is the total pseudoparticle number, zi and U{i) are the position 
and colour orientation of pseudoparticle i, respectively, and Cj = 1,1 indicates if 
pseudoparticle i is an instanton or an anti-instanton. 

To perform calculations, we take an approximate density function that describes 
the instanton liquid as an ensemble of free particles with some (unspecified) strong 

^ These values and the corresponding picture are obtained in the realistic case of three light-quark 
flavours (u, d, s): when comparing the analytic results with the numerical lattice calculations, we have 
to take into account that the latter are performed in the quenched approximation of QCD. 

^For generality, we keep A^i and Ni fixed but independent in the calculation; we will set A'l — Nj when 
needed. 
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core repulsion, which approximately localises the instantons in distinct "cells" of 
volume V^eih i-e-, 

TV 

V{{z})=M~'Y.Il^D,{^PU)), M = mvJ^,, (2.4) 

p j=l 

where {Dj} is a partition in cells of the space-time region D occupied by the 
ensemble, UjDj = D, XDj is the characteristic function of Dj^ P denotes a per- 
mutation of {1,...,A^}, and M is the normalisation; the volume of each cell is 
/ d'^zxD-j{z) = VceW- The coordinates of pseudoparticle j are denoted as Zj, and 
we take the first A'l pseudoparticles to be instantons, and the remaining Ni to be 
anti-instantons , 

= (2.5) 

[ I, A^i + 1 < yfc < iV. ^ ^ 

The function (j2.4p is by no means a precise description of the instanton ensemble: 
it has the purpose of capturing the main features of the liquid picture, namely the 
diluteness and the uniform density of the medium, while at the same time keeping 
the calculation feasible. Notice that T> does not depend on the colour orientation 
of the pseudoparticles. The precise shape of the cells should not be important, and 
one can choose it according to the geometry of the problem. 

For a local operator 0\oca\{x)-, only configurations with a single pseudoparticle 
located near the point x have both non-negligible effects and non-negligible measure 
in configuration space, due to the diluteness of the medium and to the short-range 
nature of instanton effects. The evaluation of expectation values is more complicated 
for non-local operators such as Wilson lines and Wilson loops, since they involve 
the value of the fields on a curve C. However, by properly subdividing C into small 
segments, one can exploit again the diluteness of the medium and the short-range 
nature of instanton effects, in order to consider each segment to be affected only by 
the field of a single pseudoparticle. 

As an example, consider the case of the expectation value of a rectangular Wilson 
loop W of temporal extension T and spatial width R (say, in the 1 direction). As it is 
well known, one can extract from this expectation value the static quark-antiquark 
potential Vqq using the relation [36l [37] 

(W) = ^Tr [W(^)tw(+)] = /Ce-'^^^?(-^), (2.6) 

where we have neglected the effect of the transverse "links" connecting the Wilson 
lines W^^-* , which describe the propagation of the quark and of the antiquark, and 
where /C is some proportionality constant. The effect of a single instanton located at 
z on an infinite Wilson loop, which we denote with w{z, I), is given by the following 
expression [6], 

t«(z,I)^l- J-A(z), (2.7) 
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where 

A(z) = 1 — cos a+ cos a-_ — n+ • n_ sin sin a_ , 

a t^R. ^ (2-8) 

The result does not depend on the colour orientation of the instanton, and it does 
not change if we replace the instanton with an anti-instanton. In order to calculate 
the effect of the instanton liquid, i.e., of the whole ensemble, on the Wilson loop, 
it is natural to see it as the result of a time series of interactions between the 
colour dipole and the pseudoparticles: this suggests dividing space-time into time- 

(3) 

slices of thickness 5T, and then each slice into cells of volume VceW = Vceii The 
situation is illustrated in Fig. [TJ Here 5T and V^J^ are determined by the following 
requirements: 

(3) 

1. only one pseudoparticle falls in V^J^ 5T on the average, i.e., nVceW = 1; 

(3) 

2. the spatial size V^^jj is large enough to accommodate the whole dipole, and also 
large enough so that the pseudoparticles which fall outside of the cell where 
the dipole lies will not affect it, being too distant; 

3. 6T is large enough for an (anti-)instanton to fully affect the portions of the two 
Wilson lines that fall in VceW, i-e., 5T can be considered infinite {5T S> p) as 
far as the colour rotation induced on the Wilson lines by the pseudoparticle is 
concerned, and so one can use the colour rotation angles a± given in Eq. (|2.8p . 

The region of 3-space where an (anti-)instanton can affect the operator is approxi- 
mately determined by requiring that the pseudoparticle is not too far from the two 
long sides (say, more than p), which implies 

V^2^{R + 2p){2pf- (2.9) 

in this way the two segments of the loop lying in a given slice fall in the same cell. 
With this choice we fulfill condition (2.); imposing condition (1.) we obtain 

1^ ^ ^ (210) 

p 4{R + 2p) 1 + 1' ^'-'"^ 

which is reliable for R < 2p, so that 6T/p > 5 in this regime, and condition (3.) is 
satisfied. For larger dipoles, the quark and the antiquark are more likely to interact 
with distinct instantons: the expectation value should then factorise for large R. 

At this point one can perform the average over the instanton ensemble: since the 
procedure is the same as the one adopted in [6], we will not report the calculation 
here, and we simply quote the result which will be used in the following: 

(W)iLM = e-^^-(«), V,q{R) = f d^z/\{z). (2.11) 
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This result is the same as the one found in Ref. (SU] with different, more sophisticated 
methods. As A{z) is locahsed around the position in 3-space of the quark and 
the antiquark, for large R the Wilson-loop expectation value factorises, and the 
potential becomes the sum of two constant and equal terms, which are interpreted 
as the renormalisation of the quark mass ^30j . 



3 High— energy meson— meson scattering and 
the Wilson— loop correlation functions in the 
ILM 

As has already been recalled in the Introduction, in the soft high-energy regime the 
mesonic scattering amplitudes can be reconstructed from the correlation function 
of two Euclidean Wilson loops. In this Section we give first a brief account, for the 
benefit of the reader, of the main points of the functional-integral approach to the 
problem of elastic meson-meson scattering, referring the interested reader to the 
original papers [HI [IJl [13l HU |T5] and to the book [38]. We shall use the same 
notation adopted in Ref. [28], where a more detailed presentation can be founcQ. 
We then critically repeat the calculation of the relevant correlation function in the 
ILM, and compare the result with the one found in Ref. |23| . 

The elastic scattering amplitudes of two mesons (taken for simplicity with the 
same mass m) in the soft high-energy regime can be reconstructed in two steps. One 
first evaluates the scattering amplitude of two qq colour dipoles of fixed transverse 
sizes Ri± and R2±, and fixed longitudinal momentum fractions /i and /2 of the two 
quarks in the two dipoles, respectively; the mesonic amplitudes are then obtained 
after folding the dipoles' amplitudes with the appropriate squared wave functions, 
describing the interacting mesons. The dipole-dipole amplitudes are given by the 
2-dimensional Fourier transform, with respect to the transverse distance Zj_, of the 
normalised (connected) correlation function of two rectangular Wilson loops, 

M^dd)is,t;Ri±,fi,R2±,f2) = -''^^1 d^z^e''^^-^^CMix;z±;l,2), (3.1) 

where the arguments "1" and "2" stand for "-Ri_l,/i" and "-R2±)/2" respectively, 
t = —\q±\'^ {q_L being the transferred momentum) and s = 2m^(l + coshx)- The 

correlation function Cm is defined as the limit Cm = hm Gm of the correlation 

r->oo 

■*Since no ambiguity can arise, with respect to Ref. [28^ we drop the subscript E and "tildes" from 
Euclidean quantities in order to avoid a cumbersome notation. 
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function of two loops of finite length 2T, 



yuKX^i ,z^,V,l) = — i, [6.1) 

where (. . .)m are averages in the sense of the QCD (Minkowskian) functional in- 

(T) 

tegral. Here VV| <{ are Minkowskian Wilson loops evaluated along the paths C\^2^ 
made up of the classical trajectories of the quarks and antiquarks inside the two 
mesons, and closed by straight-line paths in the transverse plane at proper times 
ibT. The partons' trajectories form a hyperbolic angle x iii the longitudinal plane, 
and they are located at (1 — fi)Ri± (quark) and —fiRi± (antiquark), i = 1, 2, in the 
transverse plane. 

The Euclidean counterpart of Eq. (|3.2p is 



GEie-, T- il; 1, 2) ^ ^'y. ''^^^f - 1, (3.3) 



where now {■■■)e is the average in the sense of the Euclidean QCD functional 

(T) 

integral. With a little abuse of notation, we denote with the same symbol W]; 2 the 
Euclidean Wilson loops calculated on the following straight-line path^, 

Ci: x1^(t) = z^ + ui^,r+(l-/i)i?i^, X^'^(r) = z^, + m^r - 

C2: X2'?(t) =^.2,.r+(l-/2)i?2M, (r) = n2^T - /2i?2M, (3.4) 

with T E [— r, T], and closed by straight-line paths in the transverse plane at 
T = ±T. Here 

til = ^sin |, 0_L, cos 1^ , U2 = ^— sin |, 0_L, cos 1^ , (3-5) 

and Ri = (0,-Ri^,0), R2 = {0, R2±,0), z = (0, zj_,0). Again, we define the correla- 
tion function with the IR cutoff T removed as C^; = lim Ge- 

It has been shown in [181 [T9| ^2] that the correlation functions in the two theories 
are connected by the analytic-continuation relation]^ 

gM{x;T;zr,l,2) = gE{-ix;iT;zr,l,2), Vx G M+, 
gEie;T;zr, 1,2) = gM{iO;-iT;zr,l,^), V0G(O,7r). (3.6) 



^The fourth Euclidean coordinate X4 is taken to be the "Euclidean time". 

^Thc functions on the right-hand side of Eqs. p.6p and p.7p are understood as the analytic extensions 
of the Euclidean and Minkowskian correlation functions, starting from the real intervals (0, tt) and M+ of 
the respective angular variables, with positive real T in both cases, into domains of the complex variables 
9 (resp. x) and T in a two-dimensional complex space. See Ref. j22| for a more detailed discussion. 
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Under certain analyticity hypotheses in the T variable, the following relations are 
obtained for the correlation functions with the IR cutoff T removed [191 122] '■ 



Cm{x;z±; 1,2) 
CE{e;zr,l,2) 



CE{-ix;z±; 1,2) 
CM{i9;z±; 1,2), 



VXGM+, 

ye e (o,7r). 



(3.7) 



We turn now to the calculation of instanton effects on the correlation function 
Ce- As it is explained in Appendix [Aj we can set /i = /2 = 1/2 without loss of 
generality; we also drop the dependence on /, in the following formulas, i.e., we un- 
derstand CE{e; z±;Ri±,R2±) = CE{e; il; /i = 1/2, i?2±, /2 = 1/2). Neglecting 
the transverse links connecting the quark and antiquark trajectories, the Wilson 
loops are written as 



>vP = ^Tr[wFtw;^)], 



i = 1,2, 



with 



Wj^^ =7'exp 



and similarly for wj'^^ Since the (long) sides of the two loops have different direc- 
tions in the longitudinal plane, their relative distance grows as we move away from 
the center of the configuration; as a consequence, only pseudoparticles falling in a 
finite interaction region have effects on both loops (see Fig. [2|), and will therefore 
contribute to the connected part of the loop-loop correlator. This region is roughly 
determined by requiring that the distance of the pseudoparticles from both loops in 
the transverse and in the longitudinal plane does not exceed ~ p. Then T^nt = ^| ^J- > 
where || and _L refer to the longitudinal and transverse planes, respectively, with 



ig / A^{X'%T))ui^dT 



i = 1,2 



(3.8) 



(3.9) 



\z±\ + 



\Ri±\ + \R2±\ 



+ 2p 



'\Ri±\ + \R2±\ 



+ 2p 



(3.10) 



(2pf 



sm ( 



(3.11) 



We now make the following approximation, considering only one instanton or anti- 
instanton in the interaction region, all the other pseudoparticles interacting with 
one loop only (or not interacting at all with the loops). Clearly, Vn blows up 
as ^ — 7- 0, vr, and so does the number of pseudoparticles in the interaction region 
"Mnt cx: np"^/ sinO; thus, the one-instanton approximation can be good if the angle 
is not too close to 0, vr. In this case we can perform the integration over the colour 
degrees of freedom independently for the two loops, except for the pseudoparticle 
in the interaction region, but since this last integration is trivial we obtain simply 



(3.12) 
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where the prime indicates that we exclude the terms corresponding to the in- 
teraction region, and refers to the interacting pseudoparticle. Here Wi{k) = 
1 — (2/A'c)Aj(rEp(^)), with Xp^^^.-^ being the position of the pseudoparticle lying in the 
k-th cell [see Eq. dlUD], and with 

Aj (x) = 1 — cos ai+ cos a^- — nj_|_ • ri.j_ sin sin aj_ , 

■^ll'^iill 



Ri 



R2 



(3.13) 



Performing the integration over the pseudoparticle positions, making use of for- 
mula (jB.7|) in Appendix [HI and dividing by the expectation values of the loops [see 
Eq. (|2.11|) ] we finally obtain for the normalised connected correlation function with 
T ^ oo 



C 



(ILM) , 



1 + n / d'^x 

'int 



Nc 



■A2(x) 



1 — n I d X —Ai{x) 



1 — n 



d^X— A2(X) 



1, (3.14) 



' mt c 

where "int" indicates that the integration range is restricted to the interaction 
region, since all the other terms cancel between the numerator and denominator. 
Expanding to first order in n, for consistency with the one-instanton approximation, 
we finally get 

)(^;z-l;i!i^,i?2±)=n(^J-) J d^x Ai{x)A2{x), (3.15) 

having extended the integration range to the whole space-time since the integrand 
is now rapidly vanishing. Exploiting the properties of the 't Hooft symbols we find 



ll"'2±||^ = 

and, moreover. 



Z_i_ ± x± 



+ I sm I - I X4 - COS I - ) xi 



_^R2± 

± 



+ I sm ( - 1 a;4 + cos ( - ) xi 



(3.16) 



71-2+ • n2- 



(£l - Xi_f 

^2 / 



R 



1± 



+ ( sm ( - 1 X4 - COS I - ) xi 



+ I sm I - I X4 + cos I - ) xi 



(3.17) 
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We can then make the dependence on 9 expHcit by performing a change of variables: 



^4 = sm ( - I X4 - cos I - ) Xl, 



x'l = sin ( - 1 X4 + cos I - ) Xl, 



(3.18) 



with 



which finally leads tc0 

cf^\0]zr^RiuR2L) = n(^^^ ^F{z^,R,^,R2±), (3.19) 

F{z^,Ri±,R2±) = y"d'xAi(x)A2(x)|e=^/2, (3.20) 

where the subscript means that the integral is evaluated setting = 7r/2 in the 
quantities (I3.16P and (I3.17p . A comparison with the result of Ref. |23j is in or- 
der. In that paper the authors perform the calculation of the correlation function 
only, but they do not take into account the effect of the pseudoparticles interacting 
with the two loops separately, which is canceled only after dividing the correlation 
function by the expectation values of the loops: their result [Eq. (63) of Ref. [23j] 
is stated in terms of an integral which is divergent, since the integrand is not van- 
ishing at infinity, while in our calculation we find a finite result. We have found 
the same analytic dependence on the angle 9 as in Ref. f23], but the prefactor 
n{2/NcfF{z±,Ri±,R2±) in Eq. 

can now be assessed numerically. 
In order to compare the analytic expression ()3.19p to the lattice data, we have to 
set the density n and the radius p to the appropriate values. The phenomenological 
estimates n = lfm~^ and p = 1/3 fm correspond to the physical (i.e., unquenched) 
case, while our numerical simulations have been performed in the quenched approx- 
imation of QCD, neglecting dynamical-fermion effects. We have then preferred to 
use the values Uq = 1.33-1.64 fm~^ for the density and pq = 0.35 fm for the av- 
erage size of the pseudoparticles, obtained directly by lattice calculations in the 
pure-gauge theory |39j . The packing fraction is a factor 1.6-2 larger than the one 
obtained with the phenomenological values, but it is still a small number (< 0.1), 
so that the dilute liquid picture is still reasonably. The ILM prediction obtained 
with these values is shown in Figs. [3H3 together with the results obtained on the 
lattice in Ref. [28j, for the loop configurations "zzz" {Ri± \\ R2± \\ z±) and "^yy" 
{Ri± II R2± -L z±), with |i2i_L| = \R2l\ = 0.1 fm. The dotted line corresponds to 



''We drop the absolute value from the Jacobian | sin 6*1, since we are limiting to 6* G (0, vr). 

^The phenomenological estimate of n is obtained assuming that instantons and anti-instantons dom- 
inate the gluon condensate G2 = {{cts / tt)F'^^,F^j^) [H [5, so that n oc 6*2. If we assumed that the same 
holds in the pure-gauge theory, we would obtain Uq/n = G2/G2 — 5.8 (using for G2 the value obtained 
in Ref. [32j, and for the quenched gluon condensate the result obtained on the lattice in Ref. |40) ) . 
i.e., a value for the pseudoparticle density Uq considerably larger than the one measured on the lattice. 
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riq = 1.33 fm~'*, while the dashed line corresponds to = 1.64 fm~^; for compari- 
son, we plot also the prediction obtained using the phenomenological values of n and 
p (solid line) . As already noticed in [28j , the functional form does not seem to be the 
correct one to properly describe the lattice data, and a second term, proportional 
to (cot 6)'^, must be added to obtain a good fit. In Tabled] we show for comparison 
the ILM prediction of the prefactor K = nq{2/Nc)'^F{z±,Ri±,R2±), and the value 
obtained with a fit to the lattice data with the fitting functions (see Ref. |28] ) 

/iLM = -f^ILM^ — 77, /iLMp = -f^^ILMp^ — 77 + -f^ILMplcot 6*)^, (3-21) 

smt^ smt^ ^ 

where /iLMp is obtained by adding the lowest-order perturbative expression [HI 
\TE\ [T9] to the ILM contribution. Notice, however, that K^i^y^p can also receive 
nonperturbative contributions from two-instanton effects |23j . 

The instanton prediction turns out to be more or less of the correct order of 
magnitude in the range of distances considered, at least around 6 = 7r/2 [where, 
as we have said before, the one-instanton approximation used to derive the result, 
Eqs. (I3.19P and (I3.20p . is expected to make sense], but it does not match the lattice 
data properly. The agreement with the data seems to be quite good at \z±\ = 0.2 fm; 
however, concerning the dependence on the relative distance between the loops, it 
seems that the ILM overestimates the correlation length which sets the scale for the 
rapid decrease of the correlation function. This is also supported by the comparison 
of the instanton-induced dipole-dipole potential, which we calculate in the next 
Section, with some preliminary numerical results on the lattice, as we show in 

Figs. EHia 



4 Wilson— loop correlation function and the di- 
pole— dipole potential 

In this Section we calculate the normalised correlation function of two (infinite) 
parallel Wilson loops, which describe the time evolution of two static colour dipoles. 
The one-instanton approximation makes no sense here, since the interaction region 
has infinite extension in the time direction, and the effect of a whole time-series 
of pseudoparticles on the two dipoles has to be considered. The above-mentioned 
correlation function can be used to extract the dipole-dipole potential Vdd by means 
of the formula |29] 

^ -TVMd,Ri,R2) Ml) 
(>Vl)(>V2) T^oo ' ^ ' ' 

where Ri and i?2 are the sizes of the two dipoles, d is the distance between their 
centers, and T is the length of the two loops. Neglecting again the transverse 
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connectors, the loops Wi^2 are written as 



with 



W 



(±) 



W, 



(±) 



■Pexp 



V exp 



Tr [W; 



-^5 



[ Ai{ut + d±Ri/2)dt 
Jo 



-ig I A4^{ut± R2/2)dt 
'o 



(4.2) 



(4.3) 



where u = (0,1), Ri = {Ri,0), d = (d, 0). To evaluate their correlation function 
[the numerator in (|4.ip ] we divide space-time into time-slices, and, in turn, the 
time-slices into cells; the procedure is the same as for the expectation value of a 

(3) 

single loop, but this time we determine the spatial size V^^^ to be large enough to 
accommodate both the dipoles. Requiring 



(4.4) 



we have to set 5T to 



6T 



4, for|d|,|i?i| ~p/2. 



6T 
P 



2.5, ioT\dl\Ri\ 



P, 



(4.5) 



which is still a fairly large value. Note however that we are overestimating the 
volume of the region of 3-space where a pseudoparticle can affect the loop, so that 
5T is actually underestimated. Again, for large distances the two dipoles interact 
with different pseudoparticles, and the correlation function is expected to factorise, 
thus giving a vanishing dipole-dipole potential as \d\ — > oo. 

Now let K be the average number of instantons that interact with the loop. This 
number grows linearly with the length T of the loop due to the uniformity of the 
liquid, and we can take, without loss of generality, T = K5T, since we are interested 
in the limit T — )• oo. We divide space-time into cells as described above, labelling 
Di, . . . , Dx the cells where the loop lives, and we define 



W^^\k) =Pexp 



W^^\k) =Pexp 



f-k6T 

/ A4{ut + d±Ri/2)dt , 

J(k-1)6T J 



-^9 



'{k-l)5T 
ckST 

-ig I Ai{ut±R2/2)dt 



{k-l)5T 



(4.6) 



and the "two-link" variables 

Ti{k) =wl-^*ik)®wl+\k), 



iij) =nT.(^), 

k=j 



(4.7) 
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which ahow us to write 

= ^ [W,(i^)],,^,^. = l-TrW,iK). (4.8) 

The integration over the colour degrees of freedom of the pseudoparticles is more 
comphcated than in the case of a single loop, although the procedure is the same. 
Starting from the K-th pseudoparticle, we have to calculate the integral 



Xi = y (iC/(/s:)TrWi(i^)TrW2(K) 



where the subscript s means that the given quantity has to be evaluated in the 
field of a pseudoparticle with standard colour orientation, where we have used the 
notation 

(4.10) 

and where 



k2 

0'^j[k[h'^k'^ + b^j[k'2^j^k[ } ^jikiSj^k2 + (^Jfc'/fe^jfe;, + "'^j[k!2^j'2k[^ ^jlk2^j2ki, 

(4.11) 



1 , 1 



(4.12) 



Performing the contractions of colour indices we obtain 

Xi = a[N^wi{K)w2{K) - w^2{K)]TtWi{K - 1)TtW2{K - 1) 

+ aN,[w^2iK) - MK)w2iK)]Wi2iK - 1) 

wher^ 

w^{k) ^ l-%{k)^,^,U = ^Tr [W;-)1"(fc)W;+)(fe)]|, = 1 - ^A,(zp(,)), 
Ai{z) = 1 — cosaj+ cosaj- — hi-^. ■ fii^ sinaj^ sinaj_, (4-13) 
"ii = ^=====, = -(d± — - z)a, = -(±— - z)a. 



^Here and in the following we adopt the following short-hand notation for quantities which depend 
on the position and the type of the pseudoparticle in the k-Va cell, f{k) — /(zp^j.), crp(j,-)). 
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and we have introduced the quantities 

Wi2(A:) = Wi{k)j,j,W2{k)j,,,. (4.14) 
The effect of a single instanton on w|^^ (k) is given by 



(+\ I exp \inf,a°'ai±\ \ , , 

U %c-2 / 



where Im is the M-dimensional unit matrix, while the effect of an anti-instanton is 
obtained by replacing rij-i- — )• —ni±. The expectation value of the quantity 

Wi2 = -^^uiK) = -^TY [W^"^%S+%S"^%^+^] (4.16) 

can be interpreted as the transition amplitude of an inelastic process, the final state 
dd being the initial one, dd, with the (say) antiquarks in the two dipoles interchanged 
(see Fig. [8]). Over the time T a certain number of such transitions can happen, and 
this is at the origin of the new terms in the formulas above. 
To proceed we also need to calculate the integral 

dU{K)Wi2{K). (4.17) 

If we write 

Wi2{K) = W,{K),,,,W2{K)j,,, = [Ti(i^)],.^,.,,.^,,j,[T2(i^)],,,.,,,.|. 



we find the same Haar integral as before, and thus, after contracting the colour 
indices, we find 



I2 = aNc[wi2{K) - wi{K)w2{K)]TiWi{K - l)Tr W2(K - 1) 

+ a[N^wi{K)w2{K) -wi2{K)]Wi2{K - 1), 

where we have introduced the quantity wi2{k), 
wi2{k) = ■^['^i{k)]j,f^j,fJs[T2{k)]j,j'^j^^^^^^ 



(4.19) 



^Tr [WS+) (A;) W;-'^ (A:) W(+) (A;) W^-^^ (A:)] U 



(4.20) 



and the quantities Wi{k), which come from the contraction of 



wi{k)w2{k) = -^[Ti(Ac)],,,.,,,.|.[T2(A;)],,,^^^^^^^ (4.21) 

c 
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and are equal to the value of the Wilson loops (of infinite length) obtained from Wj 
interchanging the antiquarks in the two dipoles, in the field of a single pseudopar- 
ticle: 



Ai (z) = 1 — cos ai+ cos a2- — ni+ • n2_ sin ai+ sin 02- 



1 - j^^iizp{k)), 



1 



(4.22) 



N, 



1 - —Aiizp^k)), 



^2(2) = 1 — cos 02+ COS ai- — n2+ • fii^ sin 02+ sinoi- 



where the functions Aj are again independent of the pseudoparticle species. Also, 
wJi2(A;) is related to Wi2{k) in the same way. 

To iterate the process we organise the previous results as follows: 



dU{K) 



TrWi(A")TrW2(A") 

Wi2(i^) 



M{K) 



TtWi{K - l)TrW2(if - 1) 
Wi2(i^-1) 



(4.23) 



where M{k) is the matrix 



M{k) 



1 



N'^wi{k)w2{k) - Wi2(fc) A^c[2i!i2(^) - wi{k)w2{k)] 



- 1 V Nc[wi2{k) - wi{k)w2{k)\ N^wi{k)w2(,k) -wuik) 
The iteration is now straightforward, and yields 



(4.24) 



NcWl2 



^j,^ TtWi{K)TvW2{K) 



'12 



1 



YlM{k) 



.k=K 



where we used the fact that 

/ TVWi(l)TrW2(l) \ 
V Wi2(l) / 



N^,WI{1)W2{1) 

Njlji2{l) 



M(l) 



Ar2 

c 
c 

Nr 



(4.25) 



(4.26) 



We are left with the integration over the pseudoparticle positions. The procedure 
is described in Appendix [Bl using Eq. (IB.7|) we then obtain 



J[M{k) 



.k=K 



1 ^ 

JJ n / (fz {viM{z, ai) + viM{z, cJj)^ 
j=K -^^i 



n 



h + n I d^z[ iyiM{z, ai) + uiM{z, ai) 



(4.27) 
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where M = M — I2 falls off to zero as \z\ — > 00, and where fj j = nj j/n. Performing 
the trivial integration over the time position, extending the spatial integration to 
the whole space, and letting (formally) K 00 with T fixed, we obtain 



where the matrix J is given by 



n ^(^) 



.k=K 



ex.p{nTJ}, 



(4.28) 



1 / N^A-B Nc{B-A) 



N^-l \ N^{B - A) N^A - B 



B 



+ 



1 / N^M-B) N,{B-A) 



^0 B ) N^-l\ Nc{B-A) N^{A-B) 
having defined 



(4.29) 



A= I (fz 

A= I d^z 

d^z 



|-)'a,(.)a.w-|-a.(.)-J-a,(,) 



(4.30) 



B -- 
B = I d^z 



Vl{Wl2{z,l) - 1) + Vl{Wi2{z,\) - 1) 
Vl{wi2{z,l) - 1) + Vl{wi2{zJ) - 1) 



To proceed further with the calculation we now set the pseudoparticle fractions to 
their phenomenological values ui = i^i = 1/2: in this case we can show that (see 
Appendix IB]) 



B = B 



2 



d-^z 



(Ai(z)A2(z) + Ai(z)A2(z) - A+(z)A_(z)) 



(Ai(z) + A2(z) + Ai(z) + A2(z) - A+(z) - A_(z)) 



(4.31) 



where the quantities 

A+(z) = 1 
A_(z) = 1 



cos Ql+ COS Q2+ 
COS Q2- COS Ql_ 



^1-1+ • '^•2+ sm ai-(- sm 02+ , 
n2- • ni_ sina2- sinai_. 



(4.32) 



are the analogues of Aj with the positions of the quark and the antiquark in the 
(say) second dipole interchanged. It is now easy to diagonalise J: denoting with 
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X = A — B, X = A — B, the eigenvalues can be written as 
A + A 1 



A± = + ^^^2 _ 1) + X ±N,^ {N^ -1){X - Xy + {X + Xry, 

_ (4.33) 
since they are differenll^, J' is diagonalisable. Denoting by II-i- the corresponding 
projectors we can write [vq = (1,0), vn^ = {N^,Nc)] 



iVc(WiW2)iLM = vo ■ {exp{nTJ}vN,) 



vo ■ 



e 



(4.34) 



and, taking the logarithm, in the large-T limit we find 

log(WiW2)iLM = nTX+ + (subleading terms). (4.35) 

In the last passage we have implicitly assumed that the projector 11+ is such that 
vo-(Il+V]\fc) > 0: we have explicitly verified that this is true in the cases that we have 
investigated numerically. Note that the correlation function of parallel loops must 
be positive, as a consequence of reflection positivity of the Euclidean theory [43]: 
this can be easily proved following Ref. |44] . Finally, recalling Eqs. (j4.ip and (|2.1ip 
we can conclude that 

Vdd{d,Ri,R2) = Vi2{d,Ri,R2) - Vgg{\Ri\) - Vg^d^l), (4.36) 

where we have set V12 = — nA+. 

In Figs. [9lll2l we show the comparison of the instanton-induced dipole-dipole 
potential with some preliminary numerical data obtained on the lattice: to our 
knowledge, lattice measurements of this quantity are not present in the literature. 
Also in this case we have performed a quenched SU{3) lattice calculation, with the 
same parameters (16^ hypercubic lattice, /3 = 6.0) used in Ref. [28], and so we 
have to use the quenched instanton density Ug and radius pq, as explained in the 
previous Section. The errors are very large for the lattice data corresponding to 
the largest distances and lengths, and a plateau has not been reached yet even at 
T = 0.8 fm at \d \ = Ofm and \d \ = 0.1 fm. However, we believe that the order of 
magnitude will not change much for larger lengths: in the various plots we show the 
available values for Vdd as obtained with loops of increasing lengths (0.4-^0.8 fm), by 
means of Eq. (j4.ip . and indeed on a logarithmic scale they seem to be quite stable. 
The lattice value of Vdd extracted from the largest-length loops is plotted at the 
correct value of the distance, while the data points corresponding to shorter lengths 
are slightly shifted, with the length increasing from left to right. The instanton- 
induced potential and the lattice data have different orders of magnitude for the 



^'^The only exception is the case X = X = 0, in which case however the matrix is already in diagonal 
form. 
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shortest distances. At larger distances there seems to be a better agreement, but the 
error bars of the lattice data are very large there, and a better precision is needed to 
quantify the importance of instanton effects at large distances. However, as already 
noticed in the previous Section, the rate of decrease with the distance seems to be 
underestimated in the ILM, which would lead to an increasing discrepancy between 
the prediction and the lattice data as the distance between the dipoles increases. 



5 Conclusions 

In this paper we have considered instanton effects, in the framework of the ILM, 
on the Wilson-loop correlation function relevant to soft high-energy scattering, and 
we have compared the results to the lattice data discussed in Ref. j28|. Using an 
approximate density function, we have critically repeated the ILM calculation of 
the relevant correlation function, already performed in Ref. [23], properly taking 
into account effects which were neglected in that paper. The analytic dependence 
on the angular variable 9 is the same as the one found in Ref. [23], which we have 
used in Ref. [28] to perform fits to the lattice data. In this paper we have instead 
performed a direct comparison of the ILM quantitative prediction, obtained by eval- 
uating numerically our analytic expression for the loop-loop correlation function in 
the relevant cases. In doing so, we have used the pseudoparticle density Uq and 
radius pq appropriate for the quenched approximation of QCD, in which the lattice 
calculation has been performed. A direct numerical comparison allows us to test 
not only the given functional form, but also the quantitative relevance of instanton 
effects on the Wilson-loop correlation function. 

The comparison of the ILM prediction with the lattice data is not satisfactory, 
although it seems to have the correct shape in the vicinity of ^ = 90°. Moreover, 
although the ILM prediction seems to be (more or less) of the correct order of 
magnitude, the rate of decrease with the transverse distance seems to be underesti- 
mated. Here one sees the importance of a quantitative comparison: terms behaving 
as 1/sin^ (which are qualitatively good in fitting the data around 6 = 90°) can 
have a different origin than instantons, and indeed it seems that the ILM is not 
able to explain them from a quantitative point of view. Also, the ILM expression 
cannot account for the small but non-zero odderon contribution to the correlation 
function which we have found in the lattice results [28j . As a final remark, note 
that, after analytic continuation into Minkowski space-time, the ILM result in the 
one-instanton approximation gives an exactly zero total cross section, since the for- 
ward meson-meson scattering amplitude turns out to be purely real. The situation 
can change after inclusion of two-instanton effects, which are expected to yield a 
constant total cross section at high energy [23]; nevertheless, an analytic expression 
for this contribution is currently unavailable. 
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In this paper we have also derived an analytic expression for the instanton- 
induccd dipole-dipole potential from the correlation function of two parallel Wilson 
loops, which we have compared to some preliminary data from the lattice. In this 
case, instanton effects seem to be negligible at short distances, less than or equal to 
the size of the dipoles, where there is a large discrepancy with the lattice results. 
The rate of decrease of the potential with the spatial distance between the dipoles 
predicted by the ILM seems also in this case to be smaller than the value found on 
the lattice; a better accuracy in the numerical calculations is needed to make this 
statement quantitative, but we think that the situation would not change from a 
qualitative point of view. 
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A Dependence of the Wilson— loop correlation 
function on the longitudinal momentum frac- 
tions 

Exploiting the symmetries of the functional integral, one ean show that the relevant 
Wilson-loop correlation function for given values of the longitudinal momentum 
fractions can always be reduced to the case /i = /2 = 1/2. We discuss here the case 
of the physical, Minkowskian correlation function Qm', the proof in the Euclidean 
case is perfectly analogous. 

As far as the dependence on the transverse vectors z_i, Ri±, and R2±, and on 
the longitudinal momentum fractions /i and /2 is concerned, the invariance of the 
theory under translations and spatial rotations implies that Qm can depend only 
on the scalar products between the relative positions of quarks and antiquarks in 
transverse space. The number of independent variables is then six, while z±, Ri± 
and fi form a set of eight variables, so two of them are redundant. Indeed, the 
loop configuration in the transverse plane is determined by the relative distances 
Ri± and R2± between quarks and antiquarks inside each meson, and by the relative 
distance between (say) quark 1 and quark 2 in the transverse plane, 

xl\T)-Xl\T) = z^ + il-fi)Ri^-{l-f2)R2±, Vr. (A.l) 
To determine which configurations are equivalent, we have to solve the equation 

^1 = ^± + {f[ - fl)Rl± (/2 - f2)R2±- (A.2) 

Since a configuration specified by {z\_, f[, f^), with z[_ given above, is equivalent to 
the configuration {z±,fi,f2) we have 

Guix; T; il; i?i±, /{, i?2±, /s) = GMix; T; zr, Ri±,h,R2±, /2); (A.3) 

one can then choose a reference value for /,', for example /( = = 1/2, and move 
all the dependence on /j inside the dependence on the impact parameter. More- 
over, since Qm enters the dipole-dipole scattering amplitude with its 2-dimensional 
Fourier transform with respect to the impact parameter, one can write 

J <fz± e"^^-^^QM{x; T; zr,R,±,h,R2±j2) 

= J dPz^ e'^^-^^QMix; T; z^ - (/i - ^)Ri± + (/s - \)R2±;Ri±, ^, i?2±, ^) 

=e^?x-[(/i-|)iiix-(/2-i)i?2x] J d^z^e'^^-^^QM{x;T;zr,Ri±,\,R2±,l). (A.4) 

It is then clear that one can consider the case /i = /2 = 1/2 only. In particular, in 
the forward case, i.e., q± = 0, the dependence on fi drops, and the integration on 
the longitudinal momentum fractions affects only the wave functions. 
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B Technical details 



B.l Integration over the position of pseudoparticles 

Consider the integral 

X= / d^'^zViiz}) n fkizp(k),'Jpik)), (B.l) 

where fk are K (possibly distinct) matrices, whose entries are functions of the 
position and of the type of the pseudoparticle lying in cell k, and Ik is the ordered 
set of the K relevant cells. We have 



1 / j4Af 



N 



P i=i 



W fk{zp{k),<^P{k)) 
k£lK 

A/'^'EK^r'' / d'^'z n XD,izpik))fkizpik),^Pik)) (B.2) 
A/'^'EKdr'' n / d^zMz,ap^,)), 



which, denoting 



becomes 



Fk{(yp(k))= I d'^z fk{z,cjp(k)), 
■Id. 



(B.3) 



where X^^^ denotes the sum over all possible sequences of K pseudoparticles, and 
r]{SK) is the number of ways in which a given sequence can be obtained from the 
ensemble, 

li and /j being the number of instantons and anti-instantons in the sequence, respec- 
tively. In the limit of large A^i, Nj and V , with the ratios vi = Ni/N, i/f = Nf/N, 
and n = N/V kept fixed, 

7]{SK)^Mn^i^l'uf, (B.6) 
and so Eq. ()B.4p simplifies to (recall that nV^eii = 1) 



2^ = ?^^E^Ni' n ^fc(^P(fc)) = '^^E n ^'^p(k)Fk{(^Pik)) 

Sk kelK Sk kelK 



n 



n E ^^^kicr) = n " E ^'^^ d*zfk{z,a). 



(B.7) 



kGlK o-=I,I 



kelK o-=I,I 
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B.2 The function B for vi = ui 

We prove that the functions B and i?, defined in Eq. (j4.30p . are equal for = z^, 
and find their exphcit form. To do so, note the fohowing about SU{2) matrices. 
One can always write u € SU (2) as 

n = u% + iu-a, n°,n^GM, {u^f + u'^ = l, (B.8) 

with = Tr2 u and u = — «Tr2(TU, with Tr2 = (l/2)Tr . Using the commutation 
and anticommutation relations satisfied by the Pauli matrices, 

[o-j,fTj] = 2ieijkcrk, {(7i,aj} = 26ijl2, (B.9) 

we can write for the product of two unitary matrices 

u = U1U2 = u^l2 + iu ■ a, (B.IO) 

with 



u'^ = [{Ulfiu2f - Ul ■ U2], U =Us+Ua, 
Us = {Uifu2 + {U2fui, Ua = -{ui A U2). 

Clearly, if we interchange ui and U2 we obtain 

V = U2U1 = v^l2 + iv ■ a, 

= U^, V = Us — Ua- 

In our case, denoting with w']j^\k) the non- trivial 2x2 part of wJ^^(A:), we have 
(suppressing the argument A;) 

wf)=«;0±).]l2 + (B.13) 



(B.ll) 



(B.12) 



so that 



where, explicitly, 



Wi = wj )^wj+) = W% + {Wi.s + Wi;a) ■ (?, (B.14) 



Wi = [(u'(_),)°(u;(+)i)° + W(^_)i ■ W(+)i], 
Wi-s = [(u'(_)i)°'W(+)j - (B.15) 

Wi-a = W{-)i A . 

Since the vectors change sign if we replace an instanton with an anti-instant- 

on, we have that under this replacement Wi-s changes sign, too, while Wi-a does not. 
Letting wi = W2 ^^w^"'''' and W2 = w^ '^^w^\ and Wj = w^'^^w- ^\ we have 

1 - 1«12 = - Tr2 [W1W2]}, 

_ 2 ^ ^ 2 (^-^6) 

l-wi2 = -^{1 - Tr2 [W1W2]} = —{1 - Tr2 [W1W2]}; 
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now, 

„o„.,o 



Tr2 [W1W2] = W1W2 - (Wi-s + Wl-a) ■ {W2-S + W2;a), ,^ ^ 

Tr2 [W1W2] = W1W2 - {Wl-s - Wi-a) ■ {W2-S " U^2;a), 



and since we are averaging with equal weights over the pseudoparticle species, we 
need to consider only terms which are symmetric under I o I, 



Tr2 [W1W2] = ivlw2 — {wi-s ■ 'W2;s + Wl-a ' 'W2;a) + antisymmetric terms 
Tr2 [W1W2] = W1W2 — {wi-s ■ W2-S + Wl-a ' W2-a) + antisymmetric terms 



Prom its definition, Eq. (j4.30p . it is now immediate to conclude that B = B. More- 
over, one can show that 

Wl:s ■ W2-S + Wl^a ' '"^2;a 

= [iw(+)lfiw(^+)2f + W(^+)i ■ W(^+-j2][iw(-)lf{w^_)2f + W^_)i ■ W^_)2] (B.19) 
- [i'W(-)lfiw(+)2f + • W;(+)2][(u'(_)2)°(w(+)i)° + W(_)2 • W(+)l], 



(B.20) 



which introducing 

w_(_ = w^'^^''^W2^^ = 'W^l2 + iw+ ■ a, 

W_ = W^ ^^W2 ^ = 'w'^12 + iw~ ■ (T, 

leads to 

Tr2 [W1W2] Isymmetric = w1w2 + W1W2 " '«^+'«^- = Tr2 [W1W2] |symmetric, (B.21) 

and thus to 

B = I dh^[{w,2{z,l) -1) + {w,2iz,l) - I)] = ^J' d'z [{A^A2 

+ A1A2 - A+A_) - (Ai + A2 + Ai + A2 - A+ - A_)] = B, 
with A± given in the text, Eq. (|4.32p . 



(B.22) 
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ij_ [fm] 


predicted 

zzz zyy 


fitted-ILM 

zzz zyy 


fitted-ILMp 

zzz zyy 


0.0 
0.1 
0.2 


0.880-1.08 0.880-1.08 
0.827-1.02 0.798-0.984 
0.692-0.853 0.607-0.748 


17.1 17.1 
7.60 5.79 
1.31 1.43 


9.86 9.86 
4.32 3.39 
0.947 1.14 



Table 1: Value of /Cilm x 10"^ for the relevant configurations: ILM prediction (first column), 
"ILM" fit (second column) and "ILMp" fit (third column). 
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Figure 1: Schematic representation of a Wilson loop in the instanton liquid. Circles 
represent pseudoparticles with the same radius p, each lying in a cell. 
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Figure 2: Longitudinal projection of the two loops, showing the partition in cells relevant 
to the case of loops at an angle 9. 
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Figure 3: Comparison of the ILM prediction, Eqs. f l3.19p - fl3.20p . with the lattice data for 
the ^""zzzj zyy^"" configuration at \z\\ = Ofm. Here and in the following figures the dotted 
line corresponds to Ug = 1.33 fm~'^ and the dashed line corresponds to Ug = 1.64 fm~^, 
with pq = 0.35 fm; the solid line corresponds to the phenomenological values n = lfm~*^ 
and p = 1/3 fm. 
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Figure 4: Same comparison as in Fig. |3]for the ^''zzz'^ configuration at \z_[_\ = 0.1 fm. 
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Figure 5: Same comparison as in Fig. [3] for the ^^zyy" configuration at \z±\ = 0.1 fm. 
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Figure 6: Same comparison as in Fig. |3]for the ^^zzz" configuration at \z±\ = 0.2 fm. 



32 



0.007 



0.006 




Figure 7: Same comparison as in Fig. [3] for the ^^zyy" configuration at \z±\ = 0.2 fm. 
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Figure 8: Schematic representation of the operators W1W2 (left) and W12 (right), cor- 
responding respectively to the processes dd — )■ dd and dd — )■ dd. White (black) circles 
represent quarks (antiquarks). 
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Figure 9: The instanton-induced dipole-dipole potential compared to the lattice data 
(on a logarithmic scale). The three lines correspond to Ug = 1.33 fm~^ (dotted line) 
and rig = 1.64 fm~'* (dashed line) with pg = 0.35 fm, and to the phenomenological values 
n = 1 fm~^ and p = 1/3 fm (solid line). Here = R2 are parallel to d, with |-Rj| = 0.1 fm. 
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